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ABSTRACT 


A computer  simulation  study  of  the  launching  of  lower  hybrid  waves 
by  external  sources  and  their  propagation  in  a finite  but  homogeneous 
plasma  slab  is  presented.  The  study  makes  use  of  a 2-1/2  dimensional, 
electrostatic  particle  code  developed  for  simulation  of  bounded  plasma. 

A number  of  antenna  configurations  are  considered:  a point  source,  two 
finite- length  capacitor  plates,  and  a phased  array  of  capacitor  plates. 
When  the  oscillating  frequency  does  not  match  a bounded  plasma  resonance, 
one  can  observe  resonance  cones,  energy  absorption  at  the  plasma  surface, 
ion  cyclotron  modulation  of  the  source,  and  energetic  ions,  depending  on 
the  parameters  chosen.  The  excitation  of  a bounded  plasma  resonance  is 
also  considered.  The  nonlinear  evolution  of  the  resonance  shows  that 
wave-particle  interactions  and  ponderomotive  force  effects  play  an 
important  role.  The  possibility  of  controlling  the  electron  temperature 
profile  is  discussed. 
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I.  INTRODUCTION 

. Plasma  heating  by  lower  hybrid  waves  has  been  considered  extensively 

1-14 

over  the  past  decade.  The  technological  advantages  of  this  scheme 

are:  (1)  RF  generators  exist  which  can  deliver  large  amounts  of  power  at 
the  frequencies  required;  (2)  the  coupling  to  fusion  devices  can  be  imple- 
mented through  waveguides  which  are  remote  from  the  neutron  environment. 
However,  if  this  scheme  is  to  be  useful,  a number  of  plasma  physics  questions 

must  be  solved.  There  is  the  problem  of  effectively  coupling  the  RF  energy 
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to  the  plasma  at  the  edge,  especially  at  the  power  levels  required.  ’ 

Further,  one  must  understand  what  happens  to  the  wave  as  it  propagates  to 

the  resonance  layer,  especially  at  large  wave  amplitudes,  where  many  non- 
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linear  effects  may  be  important.  Thirdly,  if  the  RF  energy  can  penetrate 
to  the  resonance  region,  there  is  the  problem  of  what  linear  and/or  nonlinear 

* J 8-12 

effects  occur  in  the  lower  hybrid  regions.  * Finally,  the  processes  by 

• which  the  energy  of  the  wave  is  absorbed  by  the  plasma,  and  the  subsequent 
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thermalization  of  the  energy  must  be  understood. 

Computer  simulation  can  be  a useful  tool  in  understanding  the  physics 
in  this  problem,  especially  when  nonlinear  effects  are  important  and  are 
difficult  to  treat  analytically.  Present  day  computers  are  not  large 
enough  nor  fast  enough,  however,  to  handle  all  the  different  scale  lengths 
necessary  in  simulating  the  entire  lower  hybrid  heating  scheme.  Instead, 
one  must  at  present  be  satisfied  with  simulating  certain  parts  of  the 

scheme  separately,  such  as  the  coupling  of  the  radiation  at  the  edge,  or 

. 

the  wave  conversion  in  the  plasma  interior.  One  can  nonetheless  identify 
the  important  kinetic  and  nonlinear  effects  in  each  area  and  thus  contri- 
bute to  the  understanding  of  the  overall  problem.  As  a first  step  in  a 
. systematic  study  of  the  lower  hybrid  heating  scheme,  we  present  here  a 
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computer  simulation  study  of  the  launching  and  propagation  of  lower 
hybrid  waves  by  external  sources  and  their  propagation  in  a finite  but 
homogeneous  plasma  slab  bounded  on  both  sides  by  vacuum. 

The  computer  model  and  simulation  parameters  are  discussed  in  Sec.  II. 
Non-resonant  oscillations  excited  by  various  sources  are  discussed  in 
Sec.  III.  In  Sec.  IV  the  excitation  and  nonlinear  evolution  of  a bounded 
plasma  resonance  are  studied.  Finally,  a summary  is  contained  in  Sec.  V. 


and  bounded  in  the  otheT.  In  the  non-periodic  dimension  one  can  specify 
the  values  of  the  potential  or  the  normal  electric  field  on  both  bound- 
aries. It  is  also  possible  to  specify  that  a vacuum  exists  outside  the 
plasma,  or  any  combination  of  these  conditions  on  each  boundary  separately. 
External  electrostatic  sources  (i.e.  charges  not  belonging  to  the  plasma) 
can  be  added  in  a natural  way  either  on  the  boundaries  or  in  the  plasma 
interior.  Finite-sized  (Gaussian)  particles  are  used  throughout,  and  all 
three  velocity  co-ordinates  are  retained  in  the  calculations. 

For  the  simulation  results  presented  here,  the  following  parameters 

and  conditions  were  specified:  Vacuum  boundary  conditions  were  imposed 

in  the  x direction  on  both  sides  of  the  slab.  External  sources,  described 

below,  are  specified  on  one  boundary.  The  external  magnetic  field  was 

oriented  in  the  periodic  (z)  direction  with  ratio  of  electron  cyclotron  to 

plasma  frequency  0 / u - 1,  and  ion-electron  mass  ratio  M./m  =64.  In 

e pe  1 e 

most  cases  the  ions  were  colder  than  electrons  with  Tg/T^  = 10  (except  for 
the  case  of  Fig.  4 where  Tg  = T^).  Particles  were  initially  uniformly 
distributed  within  the  slab  with  Gaussian  velocity  distributions,  and  they 
were  reflected  elastically  from  the  boundaries,  with  the  full  dynamics 
calculated  correctly  to  second  order  in  the  time  step  At  * 0.2  “pe"1* 

The  reflection  of  particles  implies  that  there  is  no  mechanism  for  net 
energy  loss  in  the  system.  The  size  of  the  system  in  the  x and  z directions 
was  Lx/ADe  = 32,  LzADe  * 128. 

The  external  source  in  these  simulations  was  determined  by  specifying 
an  external  surface  charge  density  on  one  boundary  of  the  slab,  oscillating 
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with  a specified  frequency  wq.  Because  the  magnitude  of  the  charge  density 
is  fixed  externally,  such  a source  corresponds  to  a constant  current  driver. 
In  all  cases,  the  phase  of  the  oscillator  was  chosen  so  that  the  external 
field  was  zero  at  t = 0,  and  the  amplitude  was  kept  small  so  that  the  plasma 
response  could  be  described,  at  least  initially,  by  a linear  theory.  The 
strength  of  the  oscillator  was  normalized  in  terms  of  a parameter  w,  the 
energy  of  the  source  in  vacuum  integrated  over  the  plasma  area  A divided  by 
the  electric  field  fluctuation  energy  of  the  plasma  when  in  thermal  equilib- 
rium and  free  of  external  sources, 

ELt^W2-  (1> 

' 

In  the  computer  model  for  the  parameters  chosen,  the  thermal  electric  field 
fluctuation  energy  density  is  about  100  times  smaller  than  the  particle 
kinetic  energy  density.  The  reason  for  this  particular  normalization  is 
. that  the  plasma  response  is  determined  by  the  length  and  distribution  of  the 

external  charge,  as  well  as  its  magnitude,  and  all  these  factors  are  included 
in  w. 
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III.  NON-RESONANT  OSCILLATIONS  DRIVEN  BY  AN  EXTERNAL  SOURCE 


It  is  well-known  that  the  normal  modes  of  oscillation  of  a bounded 
plasma  form  a discrete  set.16  Therefore,  the  plasma  will  respond  differ- 
ently depending  on  whether  or  not  the  frequency  and  a wavenumber  of  the 
spectrum  excited  by  the  source  match  simultaneously  the  frequency  and  wave- 
number  of  one  of  the  normal  modes  of  the  plasma.  For  perfect  matching,  the 
electric  field  for  the  single  resonant  mode  grows  secularly  with  time  until 
some  nonlinearity  stops  the  growth  . If  there  is  a mismatch,  the  plasma 
responds  with  non-resonant  driven  oscillations.  Since  a finite  size  driver 
contains  a spectrum  of  wavenumbers,  there  will  generally  be  a response  of 
many  non-resonant  modes  even  when  one  mode  satisfies  the  resonance  condition. 

The  typical  non-resonant  plasma  response  to  a number  of  "antenna  config- 
urations" of  external  sources  will  be  considered  next,  along  with  the  signifi- 
cant kinetic  and  nonlinear  effects.  Special  attention  will  be  paid  to  the 
coupling  of  the  source  to  the  plasma  at  the  surface,  the  spatial  distribution 
of  the  heating,  and  the  evolution  of  the  velocity  distribution. 

A.  Point  Source  Antenna  and  Lower  Hybrid  Resonance  Cones 

A well-known  feature  of  the  propagation  of  lower  hybrid  signals  is  the 

4 5 

existence  of  resonance  cones.  * These  cones  arise  because  for  a given 
frequency,  the  group  velocities  away  from  the  source  for  all  wavenumbers  with 
the  same  sign  are  parallel,  in  the  cold  plasma  limit.  In  order  to  clearly 
observe  resonance  cones  in  a model  which  is  periodic  in  one  direction,  it  is 
necessary  that  the  cones  from  each  of  the  periodic  sources  must  interfere 
constructively.  In  addition,  the  frequency  and  wavenumbers  of  the  source 
should  not  correspond  to  some  bounded  plasma  mode  in  order  to  avoid  a single 
mode  dominating  the  plasma  response.  Since  a resonance  cone  is  composed  of 
the  sum  of  many  Fourier  modes,  a point  source  (actually  a line  source  in 
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three  dimensions)  was  chosen  as  the  driver  in  this  case,  because  it  contains 
a broad  band  spectrum.  A driving  frequency  which  satisfies  the  above  condi- 
tions was  then  chosen  (u>o  = 2.7  w^) . The  parameter  w specifying  the 
strength  of  the  driver  was  set  equal  to  2/3. 

The  system  was  followed  for  about  11  oscillation  periods.  It  was  found 
that  the  resonance  cones  were  set  up  in  about  one  oscillation  period,  but  the 
amplitude  of  the  potential  oscillations  was  rather  weak.  In  order  to  see  the 
cones  clearly,  the  potential  at  every  point  in  space  was  Fourier  analyzed  in 
time  to  extract  the  plasma  response  at  the  driver  frequency  u)Q.  The  result  is 
shown  in  Fig.  1.  The  two  cones  are  clearly  seen,  and  make  an  angle  6 = 27° 
with  respect  to  the  magnetic  field;  however,  they  are  broader  than  expected 
from  cold  plasma  theory.  The  reason  the  cones  are  broad  is  that  the  shorter 
wavelength  components  of  the  potential  are  quickly  absorbed  by  electron  Landau 
damping  at  the  surface. 

. The  resonance  cones  shown  in  Fig.  1 can  be  Fourier  analyzed  in  the  parallel 

direction.  This  analysis  shows  that  although  the  point  source  has  a broad  spec- 
trum of  wavenumbers,  only  the  three  largest  wavelength  modes  appreciably  contri- 
bute to  the  cones  in  the  plasma  interior.  The  remaining  modes,  whose  parallel 

phase  velocities  satisfy  to  /k  < 2.5  v„.  , are  absorbed  at  the  surface.  (The 

o m - tn 

parallel  wavenumber  is  given  by  k^  = 2mn/Lz  for  integer  m and  v ^ =ATe/me'  is 
the  electron  thermal  velocity.) 

The  simulation  results  reveal  a number  of  other  important  effects.  In 
Fig.  1,  one  can  observe  the  dramatic  potential  drop  due  to  the  plasma  sheath 
surrounding  the  source.  The  effect  of  the  sheath  is  to  further  reduce  the 
penetration  of  the  source.  Examination  of  the  time  evolution  of  the  total 
electric  field  energy  wg^,  shown  in  Fig.  2 and  normalized  to  the  kinetic  energy 
. of  a thermal  electron  by 
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where  A is  the  dimensionless  Debye  length,  reveals  another  interesting 
kinetic  effect.  Superimposed  on  the  driven  oscillations,  there  is  a modula- 
tion of  the  pump  at  twice  the  ion  cyclotron  frequency  , which  manifests 
itself  as  a modulation  of  the  electric  field  energy  at  4th.  For  all  the 
non-resonant  cases  studied,  a modulation  at  the  ion  cyclotron  frequency  or 
its  first  harmonic  was  observed  whenever  the  physical  length  of  the  driver 
was  comparable  or  smaller  than  the  ion  Larmor  radius.  Although  this  modula- 
tion effect  needs  further  study,  it  is  probably  due  to  ion  kinetic  effects 
(ion  Bernstein  modes)  being  excited  in  the  sheath  region,  perhaps  due  to  the 
sudden  turning  on  of  the  pump. 

The  kinetic  energy  of  each  species  a is  normalized  in  the  same  manner 
as  the  electric  field  energy,  according  to 


De  1 r 2 
w = t-=—  -=•  m ) . v. 
a kTfi  2 a j 


(3) 


where  the  sum  is  over  the  particles  of  species  a.  The  evolution  of  the 
electron  kinetic  energy  wg  (Fig.  3)  shows  a net  energy  gain  of  about  5% 
after  11  oscillation  periods,  with  the  gain  all  due  to  the  increase  of 
parallel  velocity.  The  slight  undulations  in  the  growth  coincide  with  the 
ion  cyclotron  harmonic  modulation  noted  previously.  The  distribution  of 
energy  was  not  spatially  uniform,  however.  The  parallel  thermal  velocity 
of  electrons  confined  to  the  region  near  the  source  increased  by  20%, 
whereas  the  thermal  velocity  of  particles  on  the  opposite  side  of  the  slab 
increased  less  than  2%.  Although  a few  fast  electrons  were  formed,  the 
heating  occurred  mainly  in  the  bulk  of  the  velocity  distribution.  The  rate 
of  energy  gain  of  the  cold  ions,  also  illustrated  in  Fig.  3,  was  given  by 


dT. 

(w  T )_1  j-i 
v pe  e'  dt 
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■ = 7*10"5  . 

This  is  about  the  same  rate  observed  when  no  source  is  present.  Further,* 
the  energy  of  the  ions  stayed  spatially  uniform.  This  indicates  that  the 
energy  gain  of  the  ions  must  have  been  due  to  collisional  heating  with  the 
warm  electrons.  The  driving  frequency  in  this  case  was  too  high  for  much 
direct  ion  response. 

B.  Short  Capacitor  Plates  as  Antenna 

For  the  case  of  the  point  source,  it  is  apparent  that  electron  Landau 
damping  at  the  surface  plays  an  important  role  in  the  excitation  of  lower 
hybrid  waves  in  the  plasma  interior.  This  "filtering"  effect  of  the  plasma 
surface  is  even  more  dramatically  demonstrated  here,  where  two  short  capacitor 
plates  oscillating  out  of  phase  are  used  as  the  external  source.  The  frequency 

(u>o  = 2u)p^)  and  length  of  the  source  was  chosen  so  that  the  wavenumber  spectrum 

had  a maximum  around  those  modes  whose  parallel  phase  velocities  were  near  the 
electron  thermal  velocity  v and  the  strength  of  the  driver  was  chosen  so 
that  w = 4/3.  The  equipotential  lines  for  the  capacitor  plates  in  vacuum  are 
illustrated  in  Fig.  4(a). 

The  computer  plasma  was  followed  for  about  12  oscillation  periods,  and 
the  steady-state  plasma  response  at  the  driving  frequency  is  shown  in  Fig.  4(b). 
The  simulation  results  show  that  most  of  the  modes  are  absorbed  at  the  plasma 
surface,  and  only  the  longest  wavelength  mode  is  able  to  penetrate  to  the 
plasma  interior.  Figure  4(c),  which  shows  the  plasma  response  predicted  by 
cold  plasma  theory  for  the  excitation  of  the  single  mode  m = 1,  verifies  that 

indeed  only  that  mode  was  excited  in  the  simulation.  In  other  respects,  such 

as  the  presence  of  the  plasma  sheath  surrounding  the  exciter,  the  modulation 
of  the  source  at  the  ion  cyclotron  frequency,  and  the  evolution  of  the  plasma 
heating,  the  results  here  are  similar  to  those  of  the  point  source. 


-10- 


C.  Phased  Array  of  Capacitor  Plates  as  Antenna 

A phased  array  of  capacitor  plates  has  also  been  used  as  the  exciting 
structure.  The  frequency  (wq  = 2a)^)  and  length  of  each  element  was  chosen 
so  that  the  wavenumber  spectrum  was  sharply  peaked  around  a single  mode 
whose  parallel  phase  velocity  was  = 1.3  v^,  which  should  lead  to  a 

strong  interaction  of  the  electrons  with  the  exciting  structure.  To  enhance 
this  interaction,  the  strength  of  the  driver  was  increased  to  w = 10.  The 
equipotential  lines  for  the  phased  array  in  vacuum  is  illustrated  in  Fig.  5(a). 
The  computer  plasma  was  followed  for  12  oscillation  periods,  and  as  shown  in 
Fig.  5(b),  there  exist  no  driven  oscillations  in  the  plasma  interior. 

The  major  observations  in  this  case  are  the  familiar  electron  surface 
heating  and  the  presence  of  significant  ion  heating,  which  is  also  largest  at 
the  plasma  surface.  The  electron  kinetic  energy  wg,  illustrated  in  Fig.  6, 
shows  an  energy  gain  of  20%  after  w^t  = 300,  mostly  due  to  an  increase  in 
the  parallel  velocity.  The  energy  gain  was  localized  in  space  to  within  a 
few  electron  Larmor  diameters  of  the  source,  as  shown  in  Fig.  7. 

The  ion  heating,  on  the  other  hand,  shows  a dramatic  increase  over 
previous  levels.  The  ion  kinetic  energy  w^,  also  shown  in  Fig.  6,  increases 
by  nearly  200%  after  w^t  = 300.  Like  the  electron  case,  the  kinetic  energy 
distribution  is  inhomogeneous  in  space  and  decreases  away  from  the  source 
(Fig.  7).  However,  unlike  the  situation  for  electrons,  the  energy  gain  for 
the  ions  was  mostly  due  to  an  increase  in  velocity  in  the  x direction,  per- 
pendicular to  the  magnetic  field.  Furthermore,  as  shown  in  Fig.  8,  energetic 

ion  tails  are  formed  in  the  ion  velocity  distribution  f.(v  ).  The  reason  for 

i x 

this  difference  in  ion  behavior  is  that  the  large  perpendicular  electric 
fields  near  the  exciting  structure  (Fig.  5(b))  strongly  accelerate  the  ions 
but  do  not  affect  the  magnetized  electrons.  The  magnetic  field  would  event- 
ually return  the  ions  back  to  the  exciter  and  they  could  be  accelerated  again. 
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IV.  RESONANT  OSCILLATIONS  DRIVEN  BY  AN  EXTERNAL  SOURCE 

In  the  previous  section,  non-resonant  driven  oscillations  were  observed 
to  have  generally  low  amplitudes  inside  the  plasma.  In  this  section,  the 
enhanced  plasma  response  when  the  external  source  excites  a bounded  plasma 
resonance  is  considered.  One  example  is  studied  in  detail,  and  particular 
attention  is  paid  to  the  nonlinear  evolution  of  the  system  and  the  space  and 

|| 

velocity-space  distribution  of  the  plasma  heating. 

The  dispersion  relation  of  the  bounded  plasma  slab,  based  on  a warm 

17 

electron-cold  ion  fluid  model  was  calculated  and  is  illustrated  in  Fig.  9. 

Here  the  frequency  is  plotted  as  a function  of  the  parallel  wavenumber  and 
each  curve  corresponds  to  different  allowable  perpendicular  wavenumbers.  The 
frequency  (u)Q  = 2.3  and  length  of  two  phased  capacitor  plates  was  chosen 

so  that  the  wavenumber  spectrum  peaks  at  a normal  mode  (point  A in  Fig.  9) 
whose  parallel  phase  velocity  is  about  6 v^,  and  thus  there  would  be  little 
electron  heating  initially  by  Landau  damping.  The  source  strength  was  chosen 
to  be  w = 1. 

The  computer  plasma  was  followed  for  27  oscillation  periods.  The  time 
evolution  of  the  electric  field  energy,  illustrated  in  Fig.  10,  shows  the 
secular  resonant  growth  followed  by  absorption  of  the  wave  energy.  Then 
somewhat  later,  the  process  repeats.  The  external  source,  whose  amplitude 
is  indicated  by  the  arrow  in  Fig.  10,  was  kept  on  during  the  entire  run. 

Examination  of  the  evolution  of  the  total  kinetic  energy  of  each  species, 
illustrated  in  Fig.  11,  shows  that  the  energy  absorption  is  dominated  by 
electrons. 

A.  Electron  Dynamics 

The  rate  of  increase  of  the  electron  kinetic  energy  and  the  amplitude 
of  the  electric  field  energy  are  found  to  be  closely  correlated,  both 
reaching  their  maximum  value  at  ujpet  = 85  and  their  minimum  value  at 


-12- 


: 

j 

I 

; 

i 


I 


u>pet  = 130.  By  examining  the  electron  velocity  distribution  function,  one 
finds  that  the  increase  in  kinetic  energy  is  entirely  due  to  electrons 

A 

quickly  being  accelerated  to  large  velocities  along  the  z direction,  as  shown 
in  Fig.  12.  One  can  see  that  electrons  having  velocities  as  large  as  10  v ^ 
(that  is,  with  100  times  the  thermal  energy)  aie  created  by  w t = 220.  This 
process  continues  until  the  resonant  wave  has  lost  most  of  its  energy. 

This  behavior  can  be  understood  in  terms  of  the  concept  of  particle 
trapping.  Initially  there  were  no  electrons  moving  close  enough  to  the  phase 
velocity  of  the  wave  to  be  trapped.  As  the  potential  $o  grows  in  magnitude, 
a few  particles  can  satisfy  the  trapping  criterion 

/e? 7 
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When  this  happens,  the  wave  accelerates  the  particles  and  loses  some  energy. 
If  the  wave  damping  is  small,  the  wave  will  regain  this  energy  when  the 
trapped  particles  are  later  decelerated  by  the  wave.  However,  if  many  parti- 
cles are  suddenly  accelerated,  the  wave  may  lose  so  much  energy  that  it  can 
no  longer  trap  the  particles,  and  the  result  is  the  generation  of  fast  parti- 
cles accompanied  by  the  decay  of  the  wave.  Tail  formation  was  observed  to 
occur  in  the  simulation  when 


(6) 


indicating  that  this  acceleration  process  occurred  when  the  wave  amplitude 
was  large  enough  to  begin  to  trap  particles  in  the  bulk  of  the  electron 

velocity  distribution.  A similar  mechanism  has  been  discussed  by  Dawson  and 

e,  18 
Shanny. 
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A significant  feature  of  the  interaction  of  the  electrons  with  the 
wave  is  that  the  interaction  is  localized  in  space  because  the  eigenmode 
(not  shown)  has  large  parallel  electric  fields  near  the  boundaries.  Since 
the  electrons  are  tightly  bound  to  the  magnetic  field  lines,  it  is  not  the 
spatial  distribution  of  the  total  electric  field  for  the  eigenmode  but  rather 
the  distribution  of  the  parallel  electric  field  which  is  important  in  acceler- 
ating the  electrons,  even  though  the  magnitude  of  the  parallel  electric  field 
is  small  compared  to  the  perpendicular  electric  field.  A graph  of  the  electron 
kinetic  energy  versus  position  in  the  perpendicular  direction.  Fig.  13,  shows 
that  the  kinetic  energy  profile  in  space  follows  approximately  the  same  pattern 
as  the  parallel  electric  field  energy  for  the  eigenmode  (not  shown). 

These  results  illustrate  a possible  method  for  controlling  the  electron 
temperature  profile  in  a bounded  plasma  by  exciting  eigenmodes  whose  wave- 
lengths are  comparable  to  the  plasma  width.  Specifically,  since  a bounded 
plasma  mode  is  determined  by  the  plasma  parameters  and  geometry,  one  can  obtain 
a desired  temperature  profile  if  a mode  with  appropriate  parallel  electric 
field  distribution  is  chosen.  For  the  case  illustrated  in  Fig.  13,  about  equal 
electron  heating  on  both  sides  of  the  plasma  was  obtained,  even  though  the 
exciter  was  only  on  one  boundary.  Of  course,  only  the  parallel  energy  is 
increased. 

B.  Ion  Dynamics 

The  ions  do  not  gain  much  energy  from  the  resonant  wave,  as  one  can  see 
from  Fig.  11.  The  phase  velocity  of  the  excited  wave  is  too  fast  to  resonate 
directly  with  the  ions,  and  except  for  the  increase  of  oscillation  energy 
evident  in  Fig.  11,  the  ions  can  only  gain  energy  in  a nonlinear  fashion.  The 
key  to  understanding  the  behavior  of  the  ions  turns  out  to  be  the  ponderomot ive 
force  on  the  electrons  and  the  resulting  density  depressions.  Furthermore,  by 
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understanding  the  ion  behavior  the  reason  for  the  second  increase  in  wave 
amplitude,  shown  in  Fig.  10,  becomes  clear. 

Since  the  normal  mode  being  excited  has  a strong  spatial  distribution  of 
electric  field  energy,  the  electrons  feel  a large  ponderomotive  force.  The 
largest  component  for  electrons  is  due  to  the  parallel  gradients  of  the 
parallel  electric  field  energy.6  For  ions  the  ponderomotive  force  is  negli- 
gibly small.  The  electrons  are  therefore  pushed  away  from  regions  of  large 
parallel  electric  field,  and  the  resulting  charge  imbalance  gives  rise  to  an 
ambi-polar  potential  which  forces  electrons  and  ions  out  of  the  high  field 
regions  together,  leaving  a density  depression. 

Ions  streaming  away  from  the  high  field  regions  are  observed  in  the 

simulation  beginning  at  u t = 100,  and  there  is  a corresponding  increase  in 

pe 

ion  kinetic  energy  (Fig.  11).  The  peak  density  depression  occurs  around 

o)  t = 240,  and  a three-dimensional  view  of  the  density  (for  electrons)  is 
pe 

shown  in  Fig.  14  (averaged  over  one  oscillation  period).  The  ion  density 

looks  similar,  although  somewhat  noisier.  This  mechanism  of  ion  heating  has 

also  been  observed  in  recent  experiments  investigating  the  focus  region  of 

19 

intense  RF  radiation  above  the  lower  hybrid  frequency.  The  electrons  gain 
a negligible  amount  of  energy  in  this  manner  compared  with  the  energy  they 
gain  directly  from  the  wave. 

Because  the  system  is  periodic  in  the  z direction,  ions  from  adjacent 

high  field  regions  accumulate  between  these  regions,  until  a potential  develops 

that  opposes  further  ion  build-up.  Shortly  after,  around  u>  t = 240,  the  ion 

pe 

kinetic  energy  shows  another  jump  because  even  though  some  ions  are  slowed  by 
this  potential,  others  are  accelerated  even  more.  Thus  a slowly-varying 
potential  modulates  the  ion  heating,  and  this  shows  up  as  undulations  in  Fig.  11. 
This  type  of  ion  heating  resulted  in  fast  tails  in  the  parallel  velocity  distri- 
bution for  ions,  with  ions  having  velocities  as  large  as  5 or  6 times  the 
background  ion  thermal  velocity. 


C.  Relaxation  Oscillations 

When  the  plasma  density  relaxes  to  a nearly  homogeneous  state  (although 
fluctuations  have  increased),  around  ui^t  = 350,  the  electric  field  returns 
to  resonant  growth.  This  recurrence  can  be  described  by  the  term  relaxation 
oscillation.  Although  the  external  source  has  been  on  all  this  time,  the 
plasma  density  had  changed  so  much  that  the  source  was  no  longer  coupling  to 
a resonant  mode.  But  when  the  density  finally  relaxes,  a resonant  mode  can  be 
excited  again.  The  whole  process  we  have  been  describing  repeats  itself, 
although  the  growth  this  time  is  not  so  rapid,  since  there  are  now  particles 
traveling  fast  enough  to  interact  with  the  wave. 

Although  in  the  example  we  have  been  studying,  it  has  been  the  sudden 
acceleration  of  electrons  which  damped  the  resonant  wave,  it  is  conceivable 
that  with  different  parameters  (such  as  a resonant  mode  with  a higher  parallel 
phase  velocity),  the  density  modification  could  have  "detuned"  the  resonance 
first,  in  which  case  the  evolution  of  the  plasma  heating  could  have  been  quite 
different.  This  possibility  needs  further  study,  because  it  might  lead  to 
larger  ion  heating. 

The  plasma  state  at  the  end  of  the  computer  run  has  several  times  larger 
fluctuations  of  density  and  electric  field  than  thermal  level.  Frequency 
analysis  indicates  that  parametric  decay  to  other  modes  was  not  significant; 
for  the  two  dimensional  geometry  chosen,  however,  such  decay  was  not  expected. 
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V.  SUMMARY 


Summarizing  the  results  for  the  non-resonant  driven  oscillations  at 
frequencies  above  the  lower  hybrid,  the  plasma  response  in  the  interior  is 
relatively  small.  At  the  surface,  the  plasma  response  is  dominated  by  electrons. 
Resonance  cones  can  be  excited  in  this  model  at  certain  frequencies.  In  all 
cases  a large  potential  drop  near  the  exciter  was  observed  due  to  the  presence 
of  a plasma  sheath.  Modulation  of  the  source  at  ion  cyclotron  harmonics  by  the 
plasma  was  also  observed.  The  plasma  surface  acts  as  a filter,  absorbing  those 
wavenumbers  whose  phase  velocities  fall  in  the  electron  velocity  distribution. 
Most  of  the  plasma  heating  occurs  at  the  surface  and  takes  the  form  of  an 
increase  of  electron  velocity  parallel  to  the  magnetic  field  but  without  the 
creation  of  fast  tails  in  the  velocity  distribution.  Ions  were  appreciably 
heated  only  for  relatively  large  amplitude  exciters  by  being  accelerated  away 
from  the  source. 

The  results  for  the  resonant  driven  oscillations  can  be  described  as 
initial  growth  followed  by  relaxation  oscillations.  The  evolution  of  the 
energy  is  dominated  by  electrons.  The  resonance  is  quenched  by  sudden  tail 
formation  and  resultant  wave  damping.  Electron  heating  is  localized  in  space 
and  follows  the  pattern  of  the  parallel  electric  field  energy  for  the 
resonant  mode.  This  offers  a possibility  for  controlling  the  electron  tempera- 
ture profile  in  space  by  excitation  of  bounded  plasma  resonances.  Nonlinear 
density  changes  associated  with  the  ponderomot ive  force  are  significant.  Ions 
gain  energy  by  being  expelled  from  the  density  cavities  by  the  ambi-polar 
potential  associated  with  the  ponderomot ive  force,  as  observed  in  recent 
experiments.  When  the  density  relaxes,  the  resonant  mode  is  excited  again, 
and  the  process  repeats. 
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FIGURES 


FIG.  1 


FIG.  2 


Electrostatic  potential  associated  with  lower  hybrid  resonance 
cones  excited  by  point  source,  (ft  /w  = 1,  M./m  - 64,  T /T.  = 10, 

6 pC  16  6 1 

a>o/<»)pe  = .35,  w = 2/3,  2 x 42  x 192  particles,  32  x 128  grid, 
particle  size  ax  = a^  = 1,  XDg  = 1) 

Time  evolution  of  the  total  electric  field  energy  for  a point 
source  exciter. 


L 


FIG.  3 


FIG.  4 


FIG.  5 


Time  evolution  of  the  scaled  electron  (w„)  and  ion  (w.)  kinetic 

e l 

energies  for  a point  source  exciter. 

Equipotential  lines  for: 

(a)  short  capacitor  plates  in  vacuum; 

(b)  oscillations  driven  by  short  capacitor  plates  as 

observed  in  simulation;  (ft  / u>  = 1,  M./m  = 64, 

6 pC  1 6 

T /T.  = 1,  u>  /u  = .25,  w = 4/3,  2 x 42  x 192 
e i o pe 

particles,  32  x 128  grid,  particle  size  av  = a « 1, 

X z 

*De  “ 

(c)  oscillations  driven  by  short  capacitor  plates  as 
predicted  by  fluid  model. 

Equipotential  lines  for: 

(a)  phased  array  of  capacitor  plates  in  vacuum; 

(b)  oscillations  driven  by  phased  array  of  capacitor  plates 
as  observed  in  simulation;  (ft  /w  = 1,  M./m  = 64, 

6 pC  1 6 

T /T.  = 10,  m /w  = .25,  w = 10,  2 x 42  x 192  particles, 
el  o pe 


pe 

32  x 128  grid,  particle  size  a = a = 1,  X 


De 


1). 
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FIG.  6 


FIG.  7 


FIG.  8 


FIG.  9 


FIG.  10 


Time  evolution  of  scaled  electron  (wg)  and  ion  (vk)  kinetic 
energies  for  a phased  array  of  capacitor  plates. 

Spatial  dependence  of  the  scaled  electron  (wg)  and  ion  (w.) 
kinetic  energies  perpendicular  to  the  magnetic  field  (source 
is  at  x = 0) . 

Ion  velocity  distribution  perpendicular  to  the  magnetic  field, 
averaged  over  all  space.  Dashed  curve  corresponds  to  initial 
distribution  and  solid  curve  to  final  distribution. 

Calculated  dispersion  relation  for  a plasma  slab  with  cold  ions 
and  Lx/XDe  = 32 • 

Time  evolution  of  the  electric  field  energy  for  the  resonantly 
driven  case.  Arrow  indicates  energy  of  exciter.  (ne/u)pe  = 1, 


V\  = 64*  Te/Ti  = 10*  wo/upe 


.28,  w = 1,  2 x 42  x 192 


particles,  32  x 128  grid,  particle  size  a = a * 1,  A, 


De 


1) 


FIG.  11 


FIG.  12 


FIG.  13 


FIG.  14 


Time  evolution  of  the  scaled  electron  (wft)  and  ion  (w^)  kinetic 
energies  for  the  resonantly  driven  case. 

Electron  velocity  distribution  parallel  to  the  magnetic  field, 
averaged  over  all  space.  Dashed  curve  corresponds  to  initial 
distribution  and  solid  curve  to  distribution  after  wave  damping. 

Spatial  dependence  of  the  electron  kinetic  energy  perpendicular 
to  the  magnetic  field.  (Source  is  at  x « 0.)  Bottom  curve 
shows  energy  dependence  after  first  wave  damping  and  top  curve 
shows  energy  dependence  after  second  wave  damping. 

Electron  density  profile  at  peak  density  depression,  averaged 
over  one  oscillation  period. 


B. 

(c) 

FIGURE  4 


figure  11 
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"a  computer  simulation  study  of  the  launching  of  lower  hybrid  waves  by 
external  sources  and  their  propagation  in  a finite  but  homogeneous  plasma 
slab  is  presented.  The  study  makes  use  of  a 2-Jj  dimensional,  electrostatic 
particle  code  developed  for  simulation  of  bounded  plasma.  A number  of 
antenna  configurations  are  considered:  a point  source,  two  finite-length 
capacitor  plates,  and  a phased  array  of  capacitor  plates.  When  the  oscilla- 
ting frequency  does  not  match  a bounded  plasma  resonance,  one  can  observe  \ 
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resonance  cones,  energy  absorption  at  the  plasma  surface,  ion 
cyclotron  modulation  of  the  source,  and  energetic  ions,  depending 
on  the  parameters  chosen.  The  excitation  of  a bounded  plasma  resonance 
is  also  considered.  The  nonlinear  evolution  of  the  resonance  shows 
that  wave-particle  interactions  and  ponderomotive  force  effects  play 
an  important  role.  The  possibility  of  controlling  the  electron  temper- 
ature profile  is  discussed. 
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